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ABSTRACT 

Axisymmetric density distributions are constructed which are invisible when viewed 
from a range of inclination angles i. By adding such distributions to a model galaxy, it 
can be made either disky or boxy without in any way affecting its projected image. As 
the inclination of a galaxy decreases from edge-on to face on, the range of 'invisible' 
densities, the uncertainty in the deprojection, and the sensitivity of the deprojection 
to noise all increase. The relation between these phenomena is clarified by an analysis 
of Palmer's deprojection algorithm. 

These results imply that disk-to-bulge ratios are in principle ill-determined from 
photometry unless the disk is strong or the system is seen precisely edge-on. The 
uncertain role of third integrals in galaxies makes it unclear to what degree this inde- 
terminacy can be resolved by kinematic studies. 

Key words: Galaxies: fundamental parameters - galaxies: photometry - galaxies: 
kinematics and dynamics 



1 INTRODUCTION 

Since many galaxies are at least approximately axisymmet- 
ric and transparent, it is desirable to be able to estimate the 
three-dimensional luminosity density p(r, 9) of an axisym- 
metric galaxy from its projected surface brightness I(x,y), 
given an assumed inclination angle i between the galaxy's 
symmetry axis (the 2-axis) and the line of sight to the ob- 
server. The Richardson-Lucy algorithm (Richardson 1972, 
Lucy 1974) has been successfully used to recover distribu- 
tions p(r, 9) for a large number of galaxies with i assumed 
different from 90° (e.g., Binney Davies & Illingworth 1990; 
van der Marel 1991; Dehnen 1995). 

The uniqueness of the resulting model galaxies is un- 
clear, however. On the one hand, Rybicki (1986) gave a 
simple argument based on the 'Fourier slice theorem' that, 
for i ^ 90°, any given surface brightness distribution, I(x, y) 
must be the projection of an infinite number of luminosity 
densities p(r,9). On the other hand, Palmer (1994) proved 
that the relationship between / and p is unique for i ^ 
provided that the density p(R,9) is 'band-limited': that is 
that its expansion in Legendre polynomials Pi(cos9), 

L 

p(r,9) = J2pi(r)Pi(cos9), (1) 

1 = 

contains only a finite number of terms. Since any plausible 
luminosity density p can be approximated to sufficient ac- 
curacy by a band-limited density p, Palmer's result suggests 
that for practical purposes the relationship between / and p 
is one-to-one for all i ^ 0°. 

If two space-densities can be found which, for fixed i, 
project to the same surface brightness distribution, then the 
difference between these densities projects to zero bright- 



ness. We refer to such an 'invisible' density distribution as 
a 'konus' density. Here we demonstrate by explicit construc- 
tion of a family of konus densities that, for i ^ 90°, an infi- 
nite number of physically plausible luminosity densities are 
compatible with a given surface-brightness distribution. In 
particular, we show that one may choose between disky and 
boxy luminosity densities for the same photometric data. 

We describe how the range of possible luminosity densi- 
ties that are compatible with a given surface brightness dis- 
tribution grows continuously as the inclination varies from 
edge-on to face-on, and how this change is reflected in the 
division of the space of possible luminosity densities p(r, 9) 
into its three parts: konus densities, densities that can be 
obtained by an extension of Palmer's inversion algorithm, 
and other visible densities. 

In Section 2, we construct a family of konus densities 
and show that by adding one of these to a given depro- 
jected brightness distribution we can turn a boxy distribu- 
tion into a disky one, or vice versa. In Section 3 we show 
that Palmer's algorithm can provide a consistent deprojec- 
tion of data that do not derive from a band-limited density 
distribution provided certain conditions on the data are sat- 
isfied. These conditions will be satisfied by most data when 
the assumed inclination angle is i ~ 90°, and become more 
and more restrictive as i — > 0. 

Section 4 discusses the effects of noise in the data, the 
division of the space of possible surface brightness distribu- 
tions into ones that can be deprojected and ones that cannot, 
and the implications of the existence of konus densities for 
the detectability of low-luminosity disks in elliptical galax- 
ies. 

Section 5 sums up. In an appendix we give a group- 
theoretic proof of the result that resolves the apparent con- 
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flict between the papers of Rybicki and Palmer; namely that 
the Fourier transform of a band-limited distribution is itself 
band-limited. 



2 KONUS DENSITIES 

We now show that there exist distinct three-dimensional 
light distributions that are both physically plausible - i.e., 
are smooth and sufficiently compact - and project to the 
same surface brightness distribution. This we do by explic- 
itly calculating a family of axisymmetric luminosity distri- 
butions which are invisible when viewed at all inclination 
angles smaller than some critical value. 

We start with a review of Rybicki's (1986) application 
of the Fourier slice theorem. By writing the axisymmetric 
density p(x) in terms of its Fourier transform A(k), 



p{x) 



1 



d 3 A; A(k) exp(ifc • x) 



(2tt)3 

and projecting this along z, say, 



I(x,y) = I dz p(x) 
1 



(2*) 



(2) 



(3) 



d k A{k x , k y , 0) 8(k z ) exp(ifc • x) 



we see that the Fourier transform of the projected surface 
density I(x,y) is the two-dimensional slice A(k x ,k y ,0) of 
A(k). This simple result is known as the Fourier slice theo- 
rem. 

Consider now the projection of an axisymmetric den- 
sity distribution. From the observed image we can obtain 
the two-dimensional slice of the density's Fourier transform 
with k perpendicular to the direction of projection. How- 
ever, because of the assumed axial symmetry, all lines of 
sight which are related to the actual line of sight by a sim- 
ple rotation around the symmetry axis must give identical 
images. Therefore by symmetry the Fourier transform is 
in fact known in all parts of Fourier space that are swept 
by rotating this two-dimensional slice around the object's 
symmetry axis. 

If the density distribution is observed edge-on, this 
sweeping process covers all of Fourier space, and the three- 
dimensional density can be uniquely recovered from the pro- 
jected image. In a face-on view, only a single Fourier plane 
is known even after the sweeping process, and, as is well- 
known, the deprojection is in this case highly non-unique. 
If, however, the direction of projection is inclined with re- 
spect to the symmetry axis by an angle i ^ 90°, then the 
rotation of the Fourier slice leaves a cone around the sym- 
metry axis with half-angle (90° — i) uncovered - the so-called 
'cone of ignorance'. The surface of this cone opens to a 
plane (half-angle 90°) in the face-on case, and closes around 
the symmetry axis in an edge-on projection (half-angle 0°). 
Any density distribution derived from a Fourier transform 
which is non-zero only in the cone of ignorance will have zero 
projected brightness and thus be a konus distribution. Con- 
versely, a density distribution whose Fourier transform is 
non-zero anywhere outside the cone of ignorance will have 
non-zero projected brightness. Hence konus densities are 
precisely those densities whose Fourier transforms are non- 
zero only in the cone of ignorance. Note that a konus density 



for inclination i is also a konus density for any inclination 
i' < i, since the cone of ignorance for i lies within the cones 
of ignorance for all smaller inclinations. 

In the following we use Cartesian coordinates (x,y,z) 
and cylindrical polar coordinates (R, <f>, z) in the galaxy- 
intrinsic frame, with z along the symmetry axis, and 
(x' , y' , z') coordinates in the frame of the observer, such 
that projection is along e z i and the x = x'-axis is the line 
of nodes. The inclination angle i is measured between e z 
and e z i, the e z /-axis having direction (0, sin i, cos i) in the 
galaxy-intrinsic frame. The boundary of the cone of igno- 
rance in fc-space is given by 



\k z \/k R = \k z \j\/kl + hi = cot(90° 



\k z \ = k R tan i. 



(4) 



(5) 



The physical density corresponding to a Fourier density 
A(k) in the cone of ignorance is then 



Pk(x) 



1 



(2*) 



dk R k R 



dk z exp(i£; z 2) 



\k z | tan i 



(6) 



x / d<t>k exp[iRk R cos(<t>k - 4>)]A(k), 



where k = (k R cos <j>k, k R sm <f>k, k z ). By symmetry, A(k) 
must be independent of <j>k, so as not to generate a depen- 
dence of pk(x) on <j>. Thus the ^-integral simply evaluates 
to 2ir Jo (k R R), where Jo is the usual Bessel function. We 
will, moreover, assume that A(k) is a symmetric function 
with respect to k z = 0, 



A(k) = \ [B(k R , k z ) + B(k R , -k z )] , 
so that pk(x) becomes the cosine-Bessel transform 
1 



(7) 



dk R k R J (k R R) 



dk z cos(k z z) A(k R , k z ). 



(8) 



2.1 A specific family of konus densities 

From the discussion above it is clear that A(k) must be zero 
on and outside the cone of ignorance. So that it lead to a 
physically reasonable space density we also require it to be 
smooth, and to have certain characteristic scales. This last 
requirement derives from the fact that a konus density must 
necessarily be positive in some parts of space and negative in 
others, and after adding a konus density to a visible one we 
want the resulting density to be everywhere non-negative. 
This requires that the konus density remain finite as r — > 
and fall off sufficiently rapidly as r — > oo. Consequently, the 
konus density must have certain characteristic scales, and 
these will be inversely related to corresponding scales in its 
Fourier density A(k). 

It is likely that a general investigation of Fourier den- 
sities that are confined to the cone of ignorance will be a 
difficult numerical problem. We have therefore sought an 
analytical example and, after some experimentation with 
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Gradshteyn & Ryzhik (1980; hereafter GR), have concen- 
trated on the following family of konus Fourier densities: 



A(k) 



phere 



k 2 exp( — an) exp(— (3kR) for k > 0, 
otherwise, 



(9) 



k = \k z \ — ku tan i. (10) 

A(k) decreases smoothly to zero on the boundary of the cone 
(k = 0) and has its main contribution in a limited region of 
Fourier space around the £; z -axis, as determined by the two 
characteristic length scales a and [3. 

With this A(kn, k z ), the cosine-Bessel transform of (8) 
can be done simply. The A; z -integral gives 



dk z cos(kzz) k exp( — an) 



d 2 a cos(zku tan i) — z sin^zkn tan i) 
da 2 a 2 + z 2 

= gi(z) cos(zkR tan i) + 32(2) sm(zkR tan i), 
where we have used formula 3.893.2 of GR and 

9i{ z ) 



(11a) 



(a 2 + z 2 



(lib) 



ff2 ( z ) = t~t~, — 2^3- 

(a 2 + z 2 y j 

Note that for z ^ a the cosine and sine terms in this equa- 
tion decrease like oc z -4 and oc 2 -3 , respectively. 

Before proceeding further with the wave density in 
eq. (9), it is instructive to replace the term exp(— fiku) by 
(2A;_r) -1 5(A;_r). Then the remaining integral can immedi- 
ately be done, and we obtain the density of a plane-parallel 
sheet 



(12) 



2tt 2 (a 2 +z 2 ) 3 

The surface density obtained on integrating along any di- 
rection with inclination i' ^ 90° is 



/,(•') = / dz' Ps (z) 



dz 



7*(*) 



1 



2ir 2 cos i' 



(13) 



dz- 



(a 2 + z 2 



Only when the sheet is seen edge-on (i' = 90°) is the massive 
wire along k z part of the slice in Fourier space about which 
information is available, and only in this case will the pro- 
jected density be non-zero, with positive or negative values, 
depending on the vertical coordinate y' = z. 

We now return to the exponential wave density (9). GR, 
formula 6.751.3, give the integral 



Co 



dkn exp(— fiku) cos(zkR tan i) Jo(IcrR) 



(A + ^A 2 + B 
V2^A 2 + B 



1/2 



where 

A = R 2 - z 2 t<m 2 1 + f3 2 , 
B = 4/?Vtan 2 i. 



(14) 



(15) 



From this we may obtain the two fcn-integrals required in 
(8) as 



(16) 



Ci = dfcii &n exp(— Pkn) cos^h tan i) J (kRR) 
Jo 

dC _ 13 D 1 + ±13 z 2 tan 2 i D 2 
V2D3 

and, with a = z tan i, 



Si = I dkRkRexp(—f3kR) sm(zkRt&ni) Jo(kRR) 

dC _ ztanz(-Z)i + 4/? 2 D 2 ) 
da ~ V2D3 



(17) 



phere 



Di = A 2 + A\J A 2 + B - B, 
D 2 =A + ^y / A 2 + B, 

D 3 = (a + ^A^b) 1 ' 2 {A 2 +B^ 12 



(18) 



Thus, finally, the density corresponding to the Fourier 
density (9) may be written as 



Pk(R, z) 



1 



[g 1 (z)C 1 (R,z) + g 2 (z)S 1 (R,z)], 



(19) 



with the functions g t (z) defined by eqs (lib), and the re- 
maining auxiliary quantities specified in eqs. (14) - (18). 

This density is (i) everywhere non-singular, (ii) com- 
pact, and (iii) by construction, invisible at all inclinations 
smaller than i. To prove (i) we note from (15) that for [3 ^ 
A 2 > and B > 0, but both cannot be zero simultaneously. 
Thus the denominator D3 of both C\ and S\ is positive 
definite. Also, for a ^ the denominator of both g\ (z) 
and g 2 (z) is positive definite. To show (ii) we consider the 
asymptotic behaviour at large R and z. For R — > 00 at fixed 
z, we have A oc R 2 , B oc R° , C\ oc A~ 3/2 oc R~ 3 , Si oc R~ 3 
and, since g t are independent of R, also pk oc R~ 3 . For 
z — > 00 at fixed R, A oc z 2 , B oc z 2 , C\ oc z~ 3 , Si oc z~ 2 , 
and pk oc g 2 (z) S\(R, z) oc z~ 5 . Finally, for the special di- 
rection given by R 2 =z 2 tan 2 i, we have A = (3 2 , B = 4/3 2 R 2 , 
and with r 2 = R 2 + z 2 d oc r" 1/2 , Si oc r" 1/2 , gi oc r" 4 , 
g 2 oc r~ 3 and so pk oc r~ 7 ^ 2 . 

To confirm that this pn(R, z) is invisible in projection 
for the inclination angle i used in its construction, or indeed 
in any more face-on projection 1 with 1 < i, we have in- 
tegrated the density (19) numerically along lines of sight, 
using a Romberg midpoint rule algorithm and eq. (18) of 
Binney, Davies & Illingworth (1990). At all image points 
tried the projected density was found to be vanishingly small 
for 1 < i, and then to rise relatively rapidly from zero for 
1 > i. 

The konus densities (19) tend at large R and small z to 
constant-scale-height disks. To see this, we first note that 
for R — > 00 and fixed kR the Bessel function in eqs. (16) and 
(17) is asymptotically 

J (k R R) ~ {^—^) cos (k R R -$)+0 ( — y . (20) 



irkRR 



Thus as R — > 00, only wavenumbers near kR = (very long 
wavelengths) contribute to the konus density (19), suggest- 
ing an asymptotic connection to the plane-parallel sheet so- 
lution (12). To make this more explicit, we rewrite the konus 
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Figure 1. Contour plots in the meridional plane of two konus densities: a = f3 = 1, i = 45° (left panel); a = 0.75, f3 = 4, i = 80° (right 
panel). Contours are uniformly spaced in log( |p | ) and for negative density are dotted. 



density (19) as 



Pk(R, z) 



1 



2n 2 R 2 



dx xJo(x) exp(— flx/R) 



x ^1(2) cos ^— £ tan +(/2( 



z) sin ^"77 x tan i 



(21) 



For R — > 00, the konus density takes its largest values near 
the i?-axis (z <C R); at 2 ~ R it is smaller by a factor of order 
(ft I R) 2 . Thus we may approximate the square bracket in 
the last equation by simply gi(z), after which the remaining 
integral gives (GR, 6.623.2) 



Pk(R, 



3i(z) 
2n 2 R 2 



dx xJo(x) exp(— flx/R) 



/3 a 3 - 3az 2 
w 2 R 3 (a 2 + z 2 ) 3 



(22) 



which is the plane-parallel sheet solution (12) multiplied by a 
factor proportional to R~ 3 . Thus asymptotically for large R, 
the konus density (19) approaches a disk-like configuration 
that consists of a thin positive-density disk surrounded by a 
thicker layer of negative density. 

2.2 Some illustrations 

Fig. 1 shows contour plots in the meridional (R, z) plane 
of two members of this family of konus densities. The first 
is constructed with parameters a = fi = 1 and i = 45° in 
eq. (9), the second with a = 0.75, fi = 4, and i = 80°. Solid 
contours in Fig. 1 delineate regions in which the konus den- 
sity is positive, dotted contours signify regions of negative 
density. The structure apparent in Fig. 1 is typical of all 
the konus densities we have generated from eq. (9) - a re- 
gion of positive-density around the minor axis is followed by 
a region of negative density at intermediate latitudes, and 
then by another region of positive density around the i?-axis. 
The approach as R — > 00 of this last region to the plane- 
parallel sheet is apparent; at the largest radius plotted, the 



asymptotic density in the mid-plane, p(R,0) ~ fi j \w 2 R 3 a 3 ) 
is accurate to better than a percent. 

These density distributions project to zero surface 
brightness for all observers viewing them from inclinations 
1 < 45° and 1 < 80°, respectively - such lines of sight pass 
through regions of both positive and negative density. On 
the other hand, in an edge-on view (i = 90°) there are lines of 
sight which pass through only the region of positive density 
near the i?-axis, or only the region of negative density above 
it, depending on height z. The structures in the right-hand 
panel of Fig. 1 appear to be flattened because for high in- 
clination the cone of ignorance is narrow, allowing only for 
fairly long-wavelength components in the radial direction. 

By adding or subtracting such density components to 
a 'visible' density distribution, such as those commonly in- 
ferred for elliptical galaxies, it is clearly possible to gen- 
erate intrinsic disk-like or peanut-like components without 
altering the observed image. Ideally one would demonstrate 
this by choosing an underlying elliptical density distribution 
Pe(R, z), such as a modified Hubble profile, say, and adding 
a variety of possible konus densities. Since here we have 
fixed the functional form of the konus density (19) a pri- 
ori, we will instead vary the functional form of the elliptical 
density distribution. Fig. 2 shows three examples of such 
superpositions 



(23) 



p(R,z) = p 3 (R,z) + f p K (R,z\a,l3,i). 

In the top panel of Fig. 2 a boxy density distribution is 
obtained by subtracting a multiple of the 45°-konus density 
of Fig. 1 from the elliptical density distribution 



P3(m) 



p m 



_1 (m + ra c ) 



(24) 



Here m is the usual spheroidal radius m 2 = R 2 + z 2 / q 2 . 
This boxy distribution and the elliptical model (24) project 
to exactly the same elliptical isophotes for all inclination 
angles i' < 45°. The middle panel of Fig. 2 shows contours 
of a strongly disky intrinsic distribution, obtained by adding 
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a multiple of the a - 
halo density 



P2{m) 



p (m 2 + ml 



2, /9 = 4, 45°-konus density to the quasi- 

(25) 



Again the isophotes of this disky model are precisely ellipti- 
cal when viewed under inclination angles % < 45°. Finally, 
the bottom panel of Fig. 2 shows the superposition of the 
quasi-halo model (25) with the 80°-konus density of Fig. 1; 
this shows that invisible disklike structures can be generated 
even for high inclination % < 80°. 

It is not entirely clear which of the features of the exam- 
ples we have shown is generic and which determined by the 
peculiarities of the particular family of konus densities with 
which we have worked. Certainly, it would be straightfor- 
ward to generate konus densities with different asymptotic 
profiles from the cx R~ 3 , z~ 5 characteristic of our family: 
densities with steeper asymptotic profiles could be generated 
by either differentiating eq. (19) w.r.t. /?, or by replacing the 
exponential in eq. (9) by a Gaussian, say. Also, it is clear 
that any konus density must be made up of approximately 
conical shells of positive and negative density, and it seems 
unlikely that any has fewer than the three cones characteris- 
tic of our family. Konus densities with more cones certainly 
exist as one may demonstrate simply by adding konus densi- 
ties of the family (19) for different values of i. For example, 
Fig. 3 shows the result of adding half the 45°-density of 
Fig. I to the 80°-density of Fig. I. The combined konus 
density is invisible at i < 45°, and at r > 4 kpc it has three 
conical regions of positive density and two regions of nega- 
tive density. Adding such a density to an ellipsoidal model 
would not merely effect the transition from disky to boxy 
distributions. 

The examples of Fig. 2 are significantly influenced by 
a combination of the R~ 3 asymptotic profile and the fact 
that for R ^ fi the region of highest density around the 
i?-axis tends to a constant-scale-height disk. Between them 
these limit the amplitude / of the konus density which can 
be added to a spheroidal distribution without generating 
implausible isodensity surfaces. 

Far from the origin, the structure of a konus density will 
be dominated by the behaviour of its Fourier transform near 
the origin of fc-space. It happens that the particular wave- 
density (9) has a pronounced peak at {ku = 0, k z = 2//?), 
with the result that waves that run parallel to the 2-axis 
become dominant at large r. The tendency to a constant 
scale-height disk discussed above is a manifestation of this 
fact. In principle there is no reason why the Fourier trans- 
form of a konus density should not be zero along the A; z -axis. 
In this case the konus density would at large R be dominated 
by waves with non-zero kujk z , and would not tend to a con- 
stant scale-height disk. Some such Fourier transforms might 
lead to analytically tractable fc-space integrals, or the inte- 
grals could be done numerically. Hence it is clear that the 
range of physically reasonable konus densities is certainly 
larger than that displayed by the family (19). 

As the inclination increases towards edge-on, the angu- 
lar scales of the konus densities decrease and the form of 
the 80°-konus density of Fig. I suggests that the available 
freedom in 8p/p is less at higher than at lower inclinations. 
However, because this statement is based on the detailed 
form of our konus densities, it is not well established. We 
shall return to this point with different arguments below. 
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Figure 2. Disky and boxy systems with elliptical surface bright- 
ness contours. Each panel shows contours in the meridional plane 
for two deprojections of the same surface brightness data. The 
dotted elliptical contours correspond to one of the densities pj of 
eqs (24) or (25). The full contours show the sum of these densities 
and a konus density of plausible amplitude /. Top panel: j = 3, 
q = 0.7, m c = 2; a = f3 = 1, i = 45°, / < 0. Middle panel: j = 2, 
q = 0.7, m c = 2; a = 2, [3 = 4, i = 45°, / > 0. Bottom panel: 
j = 2, q = 0.5, m c = 3; a = 0.75, f3 = 4, i = 80°, / > 0. 
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Figure 3. A linear combination of the two konus densities of 
Fig. 1: half the left panel has been added to the right panel. The 
resulting konus density has five conical regions and is invisible at 
inclinations i < 45°. 

We have illustrated the effects of deprojection degen- 
eracy by adding konus densities to intrinsically spheroidal 
distributions, whose isophotes are elliptical. But we could 
equally well have added konus densities to disky distribu- 
tions that had pointed isophotes. Then we would have found 
that the amplitude and thickness of the disk implied by given 
data were highly ambiguous. 



3 LEGENDRE EXPANSIONS 

How do the results just presented relate to Palmer's (1994) 
proof that a band limited density distribution can be 
uniquely recovered from its projected surface brightness for 
all inclinations i ^ 0? To understand this we will first review 
Palmer's algorithm, and then return to the konus densities 
described above. We employ galaxy-intrinsic spherical polar 
coordinates (r,8,(f>). 

3.1 Palmer's algorithm 

The Legendre polynomial expansion of a typical, smooth 
luminosity density p(r, 8) will be the infinite sum 



a(r,0) = ^pi(r)Pi(cos0). 



(26) 



For simplicity we assume that p is symmetric about the 
equatorial plane, with the result that only even I need be 
included in the sum. Let (s, ip) be suitably oriented polar 
coordinates for the plane of the sky. Then the projection I 
of p may be written 



/(*,¥>) = §/„ + 



E 



I n (s) cos nip. 



(27) 



If p were band-limited, the sum over n in (27) would run only 
up to n = L, the largest value of I involved in the Legendre 



polynomial expansion of p. Thus the Fourier decomposition 
of the projection of a band-limited density distribution is 
band-limited in the conventional sense. 

We consider the band-limited approximations p ( - i "' to p 
that are obtained by setting to zero all the I n with n > L. 
Let the Legendre polynomial expansion of p ( - i "' be 



£^ L >(r)PKcos6l). 



(28) 



By expressing the (r, 8, <f>) coordinates natural to the galaxy 
in terms of coordinates (s, ip, z') natural to the observer, and 
integrating along z' , Palmer shows that 



In(s) = 4 ^ pf(cosi) / p\ L \r) Pi (p) 



r dr 



(29) 



(!-„) 



where p" is related to the conventional associated Legendre 
function P™ by 



Pi 



(l + \n\)l 



Pi" and 



cos ( arcsin(s/r)) 



(30) 



[With this definition, dx (p^) 2 = 2/(21 + 1).] Palmer 
shows that when i ^ 0, equations (29) can be converted 



into expressions for the p\ L ^ in terms of the I n , which may 
be determined observationally. Specifically, 



p\ L \r) = 
where 



2 l V.[(2l) 



H-i/2 



2%rp\ (cos %) 



d_ 

d<7 



da 



(31) 



L 

4 ^ p l y (cos i) 

l 1 = 1 + 2 



(32) 



P { v\r)Pv{^)— f 

V' 



r dr 



Palmer's proof that band-limited densities can be uniquely 
inverted follows from the fact that for i ^ this system 
of equations can be straightforwardly solved: one first finds 
P~l \ which depends only on II- Then one solves for p^ 2 > 
which depends on Il-2 and p~^\ and so on down the series 
of the p\ L \ 



3.2 Extension to non-band-limited data 

We now investigate under what circumstances Palmer's al- 
gorithmjcan be used to deproject a surface-brightness distri- 
bution I(s, (f>) that is not band-limited. Consider the differ- 
ence between the band-limited approximations p^ i+2 -' and 
p( L ) to p that are obtained with Palmer's method by trun- 
cating the Fourier series (27) for I at n = L + 2 and n = L, 
respectively. We have that 



Spi 



J L + 2 ) 



Pi -Pi 
2'n[(2Q!]- 1/2 
2%rp\ (cos %) 

where for I = L + 2 

&Il + 2 = Il + 2, 



d I81i(ar)l da 
dal a 1 \ V^2~^ 



(33) 



(34) 
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and for I < £ 

L + 2 y CO ^ 

SIi(s) = -4 y^pj,(cosi) / fipv(r)p\,(p,)- 



l' = l + 2 



(35) 



Consider now the way in which fipi depends upon the 
signal Il+2 under the assumption that the are slowly 
varying functions of s. For large I, the integral in equation 
(33) may be estimated as follows. We have 



_d_ I fili(ar) ^ ^ Ifili(ar) 
da[ a 1 \~ a^ 1 ' 



(36) 



and the dominant contribution to the integral comes from 
a ~ 1, so 



d r67i(<rr)] da 



j d<r L a 1 J V(T 2 - 1 



-iPi8I,(r] 



da 



1 a'+Wv 2 - 1 



(1-1)!! 7T - 

-' — ■r77 J — —pifih(r) 



(37) 



I!! 2' 

where /?; is a number of order unity. Thus from equation 
(33) 



\6pi\~l- 



(1-1)!! / (21)!! l3iSI,(r) 



III V ( 2l ~ !) !! 4r^(cosi')' 



(38) 



Note that since I = !!!(!- 1)!!/[(Z - l)!!(f - 2)!!)], we have for 
large I that 111/(1 — 1)11 ~ y/l. Hence for large I, 



,r / x, ,3/4 <5£(r) 

rpJ(cos z) 



(39) 



Similarly, the integral in equation (35) can be estimated as 
r dr 



fipv(r) p' v (n) 



vSpv(s) I pi'(V 1 -s 2 A 2 ) 



(40) 



Vr 2 — s 2 

, (-I)'tt / (i' + Q! (i-3)!! 
= ««- - (s) ^/ (TTTT)! (T^TT 

X3 f 2 (l^±i,_l! f i,i f i,l +1 ,I ;1 ). 

where we have used equation 7.132.6 of GR. Here 3F2 is the 
generalized hypergeometric series, which in this case termi- 
nates after i(l' — I) + 1 terms. In general we find that a 
reasonable approximation to the resulting scaling of fill is 

p t , (cos 1) s opii (s) 



\Sfi(*)\~ J2 



l' = l + 2 
L + 2 



V(l/V) 



El fl'\ p\,(cosi) CT ,. 



(41) 



i' = i + 2 



1' 1 / 4 V y pj;(cosz) 



This equation allows us to estimate the effect on the 
recovered density distribution of adding one more term Il+2 
to the expansion of the surface brightness: every coefficient 
pi depends through (31) on the corresponding 'effective data' 
coefficient £, and all of these are modified by the addition of 
Il+2 in a way that we now estimate from (41). By equation 
(34), 8Il+2 = Il+2, so from (41) we have 



SIl 



I + 2\ 2 PL + 2( cosl ) 



(£ + 2)^ V L ) p\Xl( 



(42a) 



-0.2 



-0.4 





- 


-—(50,20) 












— -^20,18) 





0.2 0.4 0.6 

x 



0.E 



1 



Figure 4. Typical Legendre functions normalized according to 
equation (30). 



Using this to eliminate SIl from the corresponding expres- 
sion for SIl— 2, we find 



SI L - 



1 



pj +2 (cosz) Pl 2 (cosi) 



(i + 2)!/ 4 i!/ 

1 



4 Pitl(cosi) 

L-2{ -\ 

p L+2 (cosi) 



(£ + 2)1/* ^(cosOJ 



p\ (cos i 

1 + 2 
1-2 



(42b) 



Il+2 



The pattern of coefficients obtained by continuing this pro- 
cess down to arbitrary I is now apparent. The dependence of 
SIi on Il+2 is governed by the ratios of the type p t ,~ 2 jp\i- 
Fig. 4 shows three typical Legendre functions. From this one 
sees that the functions p\i cx sin' i that occur in the denom- 
inators above fall monotonically as the inclination changes 
from edge-on to face-on, and become very small by a value 
of x = cos i which diminishes as I' increases. The func- 
tions p t ,~ 2 that appear in the numerators above oscillate 
for small cos i, and tend to zero ~ sin 1 -' ~ 2k ) % as cos i — ► 1. 
Consequently, the ratios of Legendre functions in (42) will 
be of order unity near i = 90° but be large for near face- 
on inclinations. In fact, in the limit i — > the products of 
Legendre functions in each term in the series of (42) will all 
be comparable because they all scale as sin -( - I ' +2-i -' i. For 
moderate inclinations and values of £, the ratio SIi/Il+2 
will be dominated by these products. Hence the addition 
of a small coefficient Il+2 can profoundly modify the effec- 
tive data £ at all orders, and hence significantly change the 
recovered density distribution. 

As we have defined it, <5£ is a function of the order £ 
of the last included term. If deprojection is to make sense, 
the difference between the deprojections that one obtains 
on truncating the data at either order £ or order £ + 2 
must tend to zero as £ — > oo. Consequently, we require 
limj,^oo SIi = for all I <C £• For i = 90°, the ratios of 
Legendre functions in equation (42) will evaluate to less than 
unity, so limj,^oo fill is zero provided the Fourier coefficients 
II fall off at least as fast as £ -2 , as they will because the 
surface brightness is a continuous function of position angle 
on the sky. However, for small i, the limit will vanish only if 
the coefficients II fall off extremely rapidly since they must 
overwhelm the very rapid growth with £ of the products of 
Legendre functions. 
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3.3 An example 

It is instructive to see how these ideas work out in prac- 
tice. Consider deprojecting a system whose noise-free sur- 
face brightness follows the modified Hubble profile 



I(x,y) 



p 2 



where 



2 



1 + 



2 i 2/2 

x +y /g 2 



(43) 



The luminosity density of this system is 



p(r,0) 



£o_ 

Hi 



where 



2 
A*3 



with 



Po 



g2-fo 

2qsro ' 



2 2 2 

q 3 = q 2 esc i 



l + ^(sin 2 + cos 2 #/gl) (44) 



(45) 



By straightforward contour integration one may show that 
for this system 

L/2 



2/n 



l 



|(g- 2 + + r 2 



yhere L is even and 

-2 



1 2 



(g 2 - 2 + l) + 2(r / S ) 2 



Since 



12 



1 + g 2 



(r > r ) 



*(1 



g 3 2 ) 



(l + ^/ g 2 sin 2 i + cos 2 
we have that II tends to zero with L as 
II 



sin 1, i(l 



2-.L/2 
93) ■ 



(46) 



(47) 



(48) 



(49) 



When this is inserted into equations (42) for SIi, the factor 
sin i+2 i in Il+2 overwhelms the factors sin -( - I ' +2-i -' i implicit 
in the products of Legendre functions, with the result that 
8Ii ~ sin' i just as /; ~ sin' i. Thus Palmer's algorithm when 
used to deproject a noise-free modified Hubble profile might 
be expected always to converge to the same density distri- 
bution, independent of the inclination at which the system 
is viewed. 

Figure 5 shows a numerical verification of this propo- 
sition. The open symbols show for a typical radius values 
of log j \pi\ that were recovered by deprojecting a modified- 
Hubble system of axis ratio qs = 0.6 and inclination i = 10°. 
The open symbols give the values obtained for L = 4, 6 and 
8, while the full symbols give the values that one obtains 
by directly expressing p(r, 8) as a Legendre series. It can 
be seen that as L is increased, good values are obtained for 
more and more of the pi. A virtually indistinguishable figure 
could be shown for any other inclination angle i > 0. Thus 
in this case, Palmer's algorithm does seem to recover the 
correct three-dimensional density from its projected surface 
brightness distribution. 



3.4 Konus densities and Palmer's algorithm 

We first demonstrate that the Legendre expansion of a konus 
density must contain an infinite number of terms; i.e., it 



-2 - 



at 
o 



-4 - 



10 



Figure 5. Legendre coefficients for a galaxy that has the modi- 
fied Hubble profile. An E4 galaxy is viewed at i = 10°. The full 
symbols show the values of log 10 \pi \ for this density distribution 
at a typical radius. The open symbols show the values of the 
same quantities that are recovered from the projected density for 
L = 4,6,8. Triangles give values of po , squares of p2, pentangles 
of pi , and so forth. 

cannot be band-limited. Indeed, if pi(r) = for 1>L, then 
by (29) it follows from the fact that Il{s) = for all s, 
that pi,(r) = for all r. Repeating this argument with 
n = L — 1, . . . in (29) we see that a band-limited density 
distribution with zero surface brightness is identically zero. 
Fig. 6 illustrates this result by showing for one of the konus 
densities of Fig. 1 the coefficients pi(r) with I < 250 at six 
values of r. 

By adding a konus density to a band-limited distribu- 
tion, we see that a band-limited surface brightness I(s,(f>) 
does not imply that the underlying density p(r, 8) is band- 
limited, while eq. (29) clearly shows that the reverse is true. 

Since, by construction, a konus density projects to zero 
surface brightness, one might think it impossible to recover 
such a density by Palmer's deprojection algorithm. Con- 
sider, however, the result of truncating the Legendre expan- 
sion of a konus density at order L. This band-limited density 
distribution projects to a band-limited surface brightness 
distribution I^ L \s,<f>). For any finite L we may in princi- 
ple deproject this by Palmer's algorithm. The band-limited 
result of this deprojection must coincide with our original 
truncated konus density, for otherwise the difference be- 
tween these two densities would be a band-limited konus 
density, which we have just shown to be impossible. Re- 
peating this experiment for larger and larger values of L, 
we obtain more and more accurate approximations to the 
original konus density, and in the limit L — > oo we obtain 
the konus density itself. Notice that this implies 



lim Pa(/' 



/ Pa( lim J 



(50) 



where Pa stands for Palmer's operator. Indeed, as L in- 
creases, I^ L ^ becomes smaller and smaller and vanishes in 
the limit L — > oo. Hence Palmer's algorithm extracts better 
and better approximations to the konus density's pi from 
smaller and smaller values of Equations (42) give some 
insight into how this is achieved, although they describe a 
limiting process that is different from the one under consid- 
eration here: now with each increment in L, every Ii varies, 
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Figure 6. The coefficients pi(r) with I < 250 for the 80°-konus 
density of Fig. 1. The larger the value of r, the more slowly p\ de- 
clines with r . This reflects the fact that at large R these distribu- 
tions approximate disks of constant scale-height. The 45°-konus 
density yields a similar figure. 

rather than just Il+2- In other words, we obtain a konus 
density from a series of Fourier series rather than from a 
single Fourier series. 

4 DISCUSSION 
4.1 Noise 

We have seen that equations (42) impose constraints on the 
behaviour of the Fourier coefficients /; at large I that must 
be satisfied if a meaningful deprojection is to be obtained by 
Palmer's algorithm. When the surface brightness distribu- 
tion I(s, (f>) is contaminated by noise, these conditions must 
be violated at sufficiently large I, since noise will tend to 
make \Ii\ approximately independent of I for large I. More- 
over, as the assumed inclination i decreases, the conditions 
imposed by equations (42) become ever more severe, with 
the result that the value of I at which a given body of noisy 
data will first violate these conditions diminishes with i. 

In practice galaxy images are deprojected using the 
Richardson-Lucy (R-L) algorithm rather than Palmer's. 
Experience shows that for a given image there is a small- 
est value of i for which the R-L algorithm yields a plau- 
sible density distribution. When using the R-L algorithm, 
both images and density distributions are usually fitted by 
low-order functions of the angles <f> and 9. This smooth- 
ing process will strongly attenuate the high-frequency power 
in both data and model, but it will not amount to strict 
truncation of the implicit Fourier and Legendre expansions. 
Consequently, the R-L algorithm will not be working with 
band-limited data or models. Equations (42) indicate that 
at sufficiently small values of i, even the residual small high- 
l terms /; in the data can make important contributions 
to the low-l coefficients pi. Since the high-order /; will be 



noise dominated, it follows that the inversion will be entirely 
noise-dominated for sufficiently small assumed inclinations. 
Clearly, the deleterious effects of noise can only be increased 
by increasing the angular resolution at which the data are 
represented. 

We have seen that the low-amplitude brightness dis- 
tributions of truncated konus densities deproject to densi- 
ties of amplitude unity. Any component of residual noise in 
smoothed data that is equal to the surface brightness of a 
truncated konus density can be matched only by projecting 
a truncated konus density. Hence no matter how such data 
are fitted, a model that fits the data accurately must include 
a truncated konus density of significant amplitude. More- 
over, the less the data have been smoothed, and therefore 
the higher the value of L at which the underlying series are 
truncated, the smaller is the amplitude of noise in the data 
that is required to produce a konus density of unit ampli- 
tude. 

It is interesting to study a practical example of a trun- 
cated konus density being introduced into a model by noise 
in the data. Fig. 7 shows a deprojection of the surface 
brightness distribution of NGC 2300 for an assumed incli- 
nation angle of i = 50°. This galaxy is approximately E2, 
and has slightly boxy (ci4 < 1%) isophotes. The deprojec- 
tion was done with the R-L algorithm of Binney, Davies 
& Illingworth (1990) as implemented by W. Dehnen. After 
10 iterations the rms deviation in surface brightness for six 
radial rays was 0.006 mag. Despite this very accurate rep- 
resentation of the nearly elliptical isophotes, the contours 
of deprojected density in the meridional plane are strongly 
non-elliptical and their shapes vary with radius. The shapes 
of these contours suggest that a truncated konus density has 
been added to a smooth elliptical model. 

Ambiguities in the result of the deprojection can be 
avoided by introducing additional constraints, such as that 
the deprojected density distribution should be as smooth 
as possible, or as elliptical as possible, etc. However, the 
realization that many elliptical galaxies display fine struc- 
ture such as weak disks casts doubts on the wisdom of such 
prescriptions. The successes of the R-L algorithm referred 
to in the introduction may in part be due to its effectively 
smoothing the angular density distribution during the de- 
projection, thus removing the higher Pi(cos9). 

4.2 What can and cannot be seen 

We have made it plausible that Palmer's algorithm can be 
applied to some non-band-limited surface brightness distri- 
butions in addition to band-limited ones. For any assumed 
inclination i, let 2Dp (i) denote the set of surface-brightness 
distributions / which yield a well-defined deprojected den- 
sity distribution in the limit L — > oo as Palmer's algorithm is 
successively applied to the Fourier decomposition truncated 
at each order L. Then 2Dp(i) comprises nearly all surface- 
brightness distributions for i = 90° and shrinks continuously 
with i until at i = it comprises only circularly-symmetric 
brightness distributions. Thus, whereas in Palmer's original 
discussion the case i = 0, in which Palmer's algorithm can- 
not be applied even to band-limited distributions, appeared 
anomalous, we now see that it is merely the endpoint of a 
continuous evolution. 

For any given value of i, there is a set 3Dp(i) of density 
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Figure 7. Contours of constant density in the meridional plane 
of a model of NGC 2300 obtained by the R— L deprojection al- 
gorithm. These are reminiscent of the sum of a konus density 
and an elliptical density distribution. The assumed inclination is 
% = 50° and the rms deviations of the projected model from the 
data are 0.006 mag. Reproduced from Hodtke (1995). 
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Figure 8. Schematic of function spaces 

distributions that can be obtained by applying Palmer's al- 
gorithm to all the members of 2Dp(i). One expects that the 
set 3Dp(i) of functions of two variables will have essentially 
the same size as the set 2Dp(i) of similar functions of which 
it is the image. No konus density lies in 3Dp(i) because to 
recover a konus density by Palmer's algorithm one has to 
change all the Fourier coefficients I n at each stage of the 
limiting process, rather than adding one more term to the 
Fourier series. Hence the set of konus densities, 3Dk(«), is 
disjoint from 3Dp(i). On adding a member of 3Dk(«) to 
a member of 3Dp (i) we obtain a density distribution that 
lies in neither set. Hence the space of all possible density 
distributions falls into three parts as is schematically illus- 
trated by Fig. 8. That figure also illustrates how 3Dk(«) 
increases from the empty set to a large part of the space as 
i diminishes from 90° to zero. 

4.3 Konus densities and disky distributions 

Konus densities are astrophysically important because ad- 
dition of a konus density can change one's model of a given 



galaxy from disky to boxy or vice versa. This strongly un- 
derlines the point made by Rix & White (1990) that many 
elliptical galaxies might contain non-negligible disks. How- 
ever, the statement of Rix & White was that these disks 
would be merely too faint to detect in currently available 
data, whereas disks contributed by konus densities could not 
be photometrically detected, even in principle. The ampli- 
tude of any disk that could be masked either by the addition 
of a konus density or by noise in the observations, decreases 
as the assumed inclination increases towards edge-on. 

One cannot place an upper limit on the luminosity of an 
'invisible' disk that can be added to a boxy galaxy without 
exploring the set of possible konus densities more thoroughly 
than we have been able to do. It seems clear that konus 
densities will always comprise a nested sequence of roughly 
conical regions of alternating positive and negative density. 
We have shown that some konus densities have more than 
the three conical regions that are characteristic of our basic 
family, but it seems unlikely that any have fewer. It is cer- 
tain that great variety is possible in the way in which the 
peak value of |p| within a conical region declines with radius 
in a konus density. For fixed inclination, our konus densities 
are characterized by two scale lengths a and fi . The latter 
determines the radius beyond which the density in the equa- 
torial plane falls as a power law, p ~ R~ 3 , while a sets the 
scale height of the disk that emerges at z <C R as R — > oo. 
These features are certainly specific to our examples rather 
than characteristic of konus densities as a whole. 

Galaxies are made up of stars moving collisionlessly on 
orbits in a given gravitational potential, and before one can 
feel free to add or subtract any component it is important 
to be assured that this can be constructed by placing stars 
on orbits. Fortunately, Lynden-Bell's (1962) demonstration 
that for any p(r, 8) there correspond infinitely many distri- 
bution functions f(E, L z ) assures us that this is so. 

In the light of the previous discussion, attempts to de- 
tect disks through their effect on kinematic observables as- 
sume greater importance. Unfortunately, uncertainty as to 
the role of a third integral in galaxy distribution functions 
makes this approach also difficult. If we knew that the dis- 
tribution function were of the simple form f(E,L z ), then 
each of the possible density distributions that are consistent 
with a given galaxy image, would predict different velocity 
profiles along the various lines of sight. Unfortunately, we 
have no reason to expect / to be of any particular form, 
and the prospects for disentangling the effects of ambiguity 
in p(r, 8) and ambiguity in phase-space structure seem slim. 
Said differently, for each of the infinitely many densities that 
are consistent with a given galaxy image, there correspond 
infinitely many possible distribution functions /, and each 
/ predicts different kinematics. It seems likely that for any 
two different densities we can find corresponding distribu- 
tion function that predict identical kinematics. 



5 CONCLUSIONS 

We have shown that for i ^ 90°, the deprojection of axisym- 
metric density distributions is non-unique in practice as well 
as in principle. That is, there are different, astrophysically 
plausible intrinsic density distributions that project to the 
same surface brightness. 
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The difference between two such densities is invisibie at 
aff sufRcientiy face-on inciinations. ft arises from a Fourier 
density in the 'cone of ignorance', the region of Fourier space 
about which the observed image contains no information. 
Exampies of such 'konus' densities can be constructed that 
are everywhere non-singuiar and decay rapidiy with radius. 

By adding a konus density to any given modei of the 
iuminosity distribution of a non-edge-on gaiaxy, one can 
change that modei from disky to boxy or vice versa without, 
even in principle, altering the fit of the model to the photo- 
metric data. At near edge-on inclination, any added konus 
density will correspond to a thin disk, but the amplitude of 
this disk could still be significant. Thus, especially in the 
case of a weak disk, the disk-to-bulge ratio is photometri- 
cally undetermined. 

Kinematic data will in some cases enable one to choose 
between luminosity distributions with identical photometric 
appearance. But it seems unlikely that even the most com- 
plete kinematic data will always resolve all ambiguity since 
the existence of a third integral for axisymmetric systems 
implies that an infinite number of distribution functions can 
generate each of the infinite number of intrinsic density dis- 
tributions that are compatible with given photometric data. 

The best prospect for making progress is to concentrate 
on studying the most edge-on systems, for which the pho- 
tometric uncertainty is least. In such systems a cold disk 
would be unmistakable but an appropriate konus density 
could add to one's model a thin hot disk without affecting 
any observable. 

We have used Palmer's deprojection algorithm to clarify 
the way in which the ambiguity of the deprojection, which 
arises from the konus densities, increases as the assumed 
inclination i decreases from edge-on. As i diminishes, the 
higher-order terms in the Fourier decomposition of the sur- 
face brightness affect more and more strongly the lower- 
order terms in the expansion of p(r, 9) in Legendre poly- 
nomials. This has two consequences. First, the range of 
surface brightness distributions that can be successfully de- 
projected diminishes as i decreases. Second, with decreasing 
i the deprojected density becomes more and more sensitive 
to noise in the data. 

Although Palmer's algorithm is helpful for understand- 
ing the deprojection problem in general terms, it is much 
harder to program and use than the R-L algorithm. Also, 
it does not constrain p to be non-negative as does the R-L 
algorithm. 

Konus densities have infinitely many terms in their Leg- 
endre expansions. The density distribution obtained by 
truncating such an expansion at high order, projects to a 
faint surface brightness distribution. When this pattern of 
surface brightness is present in smoothed noisy data, a trun- 
cated konus density of large amplitude is generated when the 
data are deprojected. We have illustrated this phenomenon 
using an image of NGC 2300. 

The results presented here for the axisymmetric case 
also have implications for general triaxial systems. For they 
strongly suggest that the range of triaxial densities that 
are consistent with given photometry is far greater than 
is commonly assumed. In particular, our results suggest 
that the three-dimensional distributions that are compatible 
with given photometry differ not only in the orientation and 
lengths of their axes, as described by Stark (1977), but also 



by the addition or subtraction of weaker, more local struc- 
tures such as disks or dumb-bells. It would be interesting to 
display explicitly triaxial analogues of konus densities. 
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APPENDIX: BAND LIMITATION OF p(K) 

Axisymmetric systems with band-limited densities p(x) 
evade the Fourier slice theorem because their Fourier trans- 
forms are also band-limited. That is, if p(k) = J d 3 x exp(ifc- 
x)p(x) is the Fourier transform of p and (k,9,(f>) are spher- 
ical polar coordinates in fc-space such that the axis 9 = 
is aligned with the galaxy's symmetry axis, then p(k,9) = 
J2i =o Pi(k)Pi(cos0). As Palmer sh ows, knowledge of band- 
limited p(k,9) in the complement of Rybicki's "cone of ig- 
norance" suffices to determine p(k,9) for all k. In fact, p(k) 
is determined by its value at only a finite number of values 
of 9 at each value of k. 

Palmer effectively demonstrates that p is band-limited 
by explicit calculation. Here we show that this result follows 
easily from a general group-theoretic argument. 

Let the functions f m (x) for m = —I, . . . , I form a basis 
for the spin-l irreducible representation of SO(3), where the 
action on /; of a rotation R £ SO(3) is 

R:f^f with f(x) = f(R T ■ x). (Al) 

(Here R is the rotation matrix associated with R ) Then by 
the uniqueness of the spin-l representation, any f m can be 
written as 

i 

f m (x)= fmm'(r)Yr(9A)- (A2) 

m' = -l 
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Fourier transforming equation (Al) we have 
f'(k) = j d 3 x' exp(ifc • x')f'(x') 




d 3 x'exp(ifc • x')f(R T ■ x') 
d 3 x exp[ifc • (R ■ x)]f(x) 
d 3 k exp[i(i? T • k) ■ x]f(x) 



(A3) 



= f(R T -k). 



Consequently, the Fourier transforms f m also form a basis 
for a (2l+l)-dimensional irreducible representation of SO(3), 
and by the uniqueness of the spin-l representation we have 



A band-limited p(x) is a function that can be written as a 
linear combination of functions that are members of spin-l 
representations of SO(3) for I < L. By the above it immedi- 
ately follows that p(k) has this last property and so is itself 
band-limited. 



This paper has been produced using the Blackwell Scientific 
Publications T^X macros. 
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